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Abstract 

Excitation-energy distributions have been derived from measurements of 5.0-14.6 
GeV/c antiproton, proton and pion reactions with ^^^Au target nuclei, using the 
ISiS 47r detector array. The maximum probabihty for producing high excitation- 
energy events is found for the 8 GeV/c antiproton beam relative to other hadrons, 
^He and p beams from LEAR. For protons and pions, the excitation-energy distri- 
butions are nearly independent of hadron type and beam momentum above about 
8 GeV/c. The excitation energy enhancement for p beams and the saturation ef- 
fect are qualitatively consistent with intranuclear cascade code predictions. For all 
systems studied, maximum cluster sizes are observed for residues with E*/A ~ 6 
MeV. 
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Much effort in nuclear pfiysics fias been devoted to tlie study of tlie formation 
of multiple complex fragments (3< Z <16), or multifragmentation, and its 
possible link to the nuclear liquid-gas phase transition [0. In particular, the 
strong interest created by the measurement of a latent heat by Pochodzalla 
et al. PI , which would signal a first order phase transition. However, in many 
cases multifragmentation of heavy nuclei is not only driven by the its thermal 
and Coulomb properties, but also by collective (dynamical) properties of the 
chaotic systems formed in energetic projectile-target interactions. 

The thermal features of the breakup process are isolated most transparently 
in reactions induced by hadron and light-ion beams at energies in excess of 2 
GeV0-§. Transport calculations predict that such beams heat nuclei rapidly, 
r < 30-40 fm/c, at the same time producing little compression, low average 
angular momentum and minimal shape distortions [^,0. While excitation en- 
ergy deposition E* in GeV hadron- induced reactions is significantly less than 
the total available energy[^H, transport calculations indicate that residues 
can be created in such collisions with E* values well in excess of the multi- 
fragmentation threshold, E* j A ~ 5 MeV||^-[l0[|. The objective of the present 



study is to investigate the relative effectiveness of various hadron beams and 
momenta for producing E* values in excess of ~ 5 MeV/A, thereby identifying 
the optimum system for studies of thermal multifragmentation and underlying 
phenomena e.g. the liquid-gas phase transition. 



In Fig. 0, excitation-energy predictions of the Toneev|TT| intranuclear cascade 
calculation (INC) are shown. Here the average excitation energy imparted to 
a ^^^Au nucleus by proton, vr" and antiproton beams is plotted as a func- 
tion of beam momentum. For p and 7r~ beams, the (£^*) values are predicted 
to be nearly identical; above about 8 GeV/c there is little dependence on 
beam momentum. These features have been verified qualitatively in charged- 
particle multiplicity studies by Hsi et al. |12| and in earlier inclusive studies by 
Porile[|13|. However, the antiproton predictions exhibit a significant increase 
in average excitation energy. This enhancement derives from reabsorption of 
some fraction of the annihilation pions (< n^ >~ 5), which complements 
the internal heating caused by the cascade of hadron-hadron collisions and 
A resonance excitations P JTT|JT^ . From the point of view of multifragmenta- 
tion, this greatly enhances the probability for forming residues excited above 
the multifragmentation threshold, as shown in the inset of Fig. 1, where E* 
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distributions for 8 GeV/c tt and p are compared. 

In this letter we present results for excitation-energy distributions derived 
from bombardments of ^^^Au nuclei with 5.0-14.6 GeV/c hadron beams. Ex- 
clusive charged-particle multiphcities and energy spectra were measured at the 
Brookhaven AGS accelerator with the Indiana Silicon Sphere, a Air detector 
array with 162 gas-ion-chamber/silicon/CsI telescopes]^. Two experiments 



were performed. The first used untagged negative beams (largely tt") at 5.0, 
8.2, and 9.2 GeV/c and positive beams (primarily protons) at 6.2, 10.2, 12.8 
and 14.6 GeV/c. The second was performed with a tagged 8.0 GeV/c nega- 
tive beam, which provided simultaneous measurement of the vr" (98%) and p 
(1%) reactions. The results at 8 GeV/c for the two vr" experiments are found 
to be identical within error bars. Further experimental details can be found 
in [|l2| , p!5| , |l6| . As part of our analysis, we have also compared with data from 
the 4.8 GeV ^He + ^^^Au reaction § at LNS Saclay and the 1.2 GeV/c p + 
^^^Au reaction from LEAR[|]. 

In the experiments reported here, calculation of the excited residue charge 
and mass was made by subtracting fast cascade particles from the target 
charge and mass using the same procedure as in ref @]. Excitation-energy 
reconstruction was performed for each event by calorimetry according to the 
following prescription: 

Mc 

E* = Y,K, + M^<K^>+Q + E^. (1) 

Here, Ki is the kinetic energy of each charged particle in an event of multi- 
plicity Mc, M„ and {K^) are the multiplicity and average kinetic energy of 
neutrons, Q is the mass difference of the reconstructed event, and E^ is a small 
term to account for gamma de-excitation of the residual nucleus and excited 
fragments. This procedure is similar to those employed in P,|l6 i p!7 i Pi[T8|-pO[] and 



a full paper on the dependence of E* values on the various assumptions of the 
reconstruction and the effects of fluctuations is in preparation. 

For the present measurements the two most sensitive parameters of the recon- 
struction procedure involve the assumption concerning {Kn) and the definition 
of the energy range for thermal-like particles. Since neutrons are not measured 
in ISiS, we use the neutron-charged particle correlations reported for LEAR 
data by Goldenbaum et a/. P]. This correlation is in good agreement with sim- 
ilar results from heavy-ion reactions[^ and model simulations P,P,p2| . How- 



ever, systematic uncertainties arise when the kinetic energies are assigned to 
the neutrons. The neutron average kinetic energies were taken from the pre- 
dicted correlation between {Kn) and E*/A by SMMp. Then Eq. 1 is iterated 
to obtain self-consistency. Similar values are obtained from an iterative proce- 
dure using E* = aT"^ and an initial value of (Kn) = T. Comparison between 
unfiltered and filtered simulations with SMM[R[ and the evaporation code SI- 



MON H] both show that the use of 3T/2 (initial {Kn) = 3T/2) overpredicts 
the values of (Kn) by as much as 30% at high E*, resulting in an overestima- 
tion of E* of about 10-12% (~ 20% for initial {K^)=2T). 

Thermal-like charged particles are defined by the spectral shapes, from which 
an upper cutoff of 30 MeV for H and 9Z + 40 MeV for heavier fragments was 
assigned [Q^. Our definition oi E* is a conservative one; for example, with the 
expanded charged particle acceptance {E/A < 30 MeV) of Hauger et a/.|]2D|, 
we obtain E* values about 25% higher than reported here. 

Finally, the stability of the E* reconstruction procedure has been tested us- 
ing SMM and SIMON. Both models give a strong linear correlation between 
the unfiltered and filtered E*. The average values are recovered by the method 
with deviation not larger than ±10% over the useful range of the data {E*/A=2- 
9 MeV). A detailed comparison will be presented in a long paper. However, 
the approach taken above should be viewed as a conservative one, as should 
the E* distributions. 

In Figure 2 we show the reconstructed probability distributions for excitation 
energy and residue mass for the range of systems studied in this work, where 
J2i P{El) = 1. Values of E* below 250 MeV become highly uncertain due 
to the dominance of neutron emission (unmeasured in ISiS) at low excitation 
energies and the suppression of M^ < 2 events by the ISiS trigger. Data for 6.2 
GeV/c and 12.8 GeV/c protons (not shown) are similar to other proton and 
pion data in Figs. 2-5. Fig.^ demonstrates that the largest population of high 
excitation energy events is achieved with the 8.0 GeV/c p beam and the lowest 
with the 5.0 GeV/c vr" beam. The 12.8 GeV/c proton distribution is slightly 
higher than for 14.6 GeV/c p, while that for the 6.2 GeV/c protons is slightly 
lower than the 8 GeV/c tt^. Thus, the data and the INC predictions of Fig. 
|I| are in qualititive agreement. Quantitatively, however, the INC calculations 
predict E* distributions that extend significantly beyond the data, as discussed 
in^ 



The residue mass distributions show a somewhat different pattern. In this case 
the 14.6 GeV/c proton beam produces the lightest residues and the 5.0 GeV/c 
71" the heaviest. This mass dependence can be understood as a consequence of 
the fast cascade, which produces an increasing number of fast knockout par- 
ticles as the beam momentum increases, and also leads to a density-depleted 
residue [|],0. This process produces the saturation in average excitation energy 
shown in the INC calculations of Fig.|I| and the data in Fig.^ That is, the in- 
crease in total available beam energy for E* deposition is counterbalanced by 
loss of energy due to mass loss AA during the fast cascade. This mass loss, 
derived from the data, is shown in the top panel of Fig.^ as a function of 
deposited excitation energy. 



The relative effectiveness of various beams in depositing high excitation en- 
ergies is shown in the bottom panel of Fig. ^ and summarized in Table 1. 
Included here are comparable data from the 4.8 GeV ^He + ^^'^Au reaction^ 
and the 1.2 GeV p + ^^^Au reaction [0]. In order to emphasize the probabil- 
ity for forming highly excited systems, we examine the ratio of total events 
with E* greater than a given value to that for events with E* > 400 MeV. 
At ii^*=400 MeV, the event reconstruction should provide the greatest self- 
consistency among the data sets. 

Figure ^ and Table |1| confirm that the 8.0 GeV/c antiproton beam produces 
a significant enhancement of high excitation energy events, particularly in 
the multifragmentation regime above 800-1000 MeV. In Table |I], the yield 
of events with excitation energy above the multifragmentation threshold for 
Au-like nuclei (a range that spans 800-1000 MeV or about 5 MeV/nucleon) 
is listed, compared to total events above E* >400 MeV {E*/A >2 MeV). 
The enhancement for the 8 GeV/c p beam is approximately 25% greater than 
the next most effective beam, 12.8 GeV/c protons. Relative to the p studies 
with 2.1 GeV/c p at LEAR |^, where negligible multifragmentation yield was 
observed, the probabilities for high E* events at 8 GeV/c are over an order of 
magnitude greater. In this regard, we note that for the 1.8 GeV ^He + ^^^Au 



system |23], for which the charged-particle multiplicity data are very similar 



to those with LEAR beams |2^, the total cross section for events with four or 



more Z > 3 fragments is only 3.5 mb. At the higher energy of 4.8 GeV, this 
cross section has grown to 83 mb. Thus, when account is made for the rapid 
growth of multifragmentation cross section with increasing beam momentum, 
the ISiS results and those of Ref. appear to be self-consistent. 

In Fig. ^ the excitation energy distributions are plotted as a function of E*/A 
of the residue. We note here, as well as in Fig ^, that events with E*/A > 
9 MeV comprise less than 1% of the data set. The same general features 
persist as in Fig. ^, except in this case the 14.6 GeV/c proton beam yields 
comparable probabilities in the region beyond about E* / A > 9 MeV. Two 
factors account for this. First, the average residue mass is lighter for reactions 
at this momentum, as shown in Fig. ^ Second, the number of events obtained 
with the 8.0 GeV/c p beam (~ 25,000) were about two orders of magnitude 
lower than for the other beams, creating larger statistical uncertainties at the 
extremes. 

Finally, in Fig. |^ we examine the dependence of the fragment size distribu- 
tions on E*/A, of relevance to discussions of critical phenomena and a nuclear 
liquid-gas phase transition [^,^]. The top panel shows the number of ob- 
served IMFs per residue nucleon and the corresponding filter-corrected value 
as a function of E*/A for the various reactions. This ratio is nearly identical 
for all systems, increasing systematically with increasing E*/A up to 9 MeV 
and then becoming roughly constant thereafter. The same uniformity is ob- 



served in the fragment charge distributions, shown in the lower panel of Fig. 
^, where the parameter r from power-law fits to the charge distributions are 
plotted as a function of E* jA. Values of r decrease steadily as the system is 
heated i.e. the probability for forming larger fragments increases. A minimum 
is reached at r ~ 2 near E* /A r^ Q MeV, followed by a slight increase (smaller 
fragments). This signifies that maximum cluster sizes are obtained very near 
the multifragmentation threshold. Thereafter, additional excitation appears to 
produce a hotter environment, leading to an increased yield of lighter particles 
and clusters. 

In summary, the heat content {E*) of equilibrium-like heavy residues formed 
in 5-15 GeV/c hadron-induced reactions has been investigated. The antipro- 
ton beam is found to be most effective in creating highly excited residues, in 
qualitative agreement with INC predictions. Relative to the threshold for mul- 
tifragmentation in such systems {E* ~ 800-1000 MeV), the enhancement of 
high excitation energy events with antiprotons is at least 25-35% greater than 
other hadrons and over an order of magnitude greater than antiprotons from 
LEAR. Above momenta of about 8 GeV/c the probability for E* deposition 
with hadron beams is nearly independent of hadron type or beam momen- 
tum, again consistent with INC calculations. The observed average number 
of IMF per residue nucleon and the power-law fits to the charge distributions 
show a universal behavior as a function of E* /A. This independence of the 
final multifragmentation state on collision dynamics suggests equilibrium-like 
behavior in the breakup of the hot residues I^S],^. Therefore hadron, espe- 
cially antiprotons around 6-8 GeV/c and pions or protons above 10-12 GeV/c, 
are very well suited to study thermal multifragmentation. For a given beam 
momentum, they provide a wide range of thermal energy {E* /A), which is an 
essential quantity for investigating latent heat in nuclear matter and related 
properties. 
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Fig. 1. Intranuclear cascade predictions |11] of the average excitation energy for 
events with E* > 50 MeV are shown as a function of momentum for p, tt~ and p 
beams incident on ^^^Au. Inset compares the excitation energy probability distri- 
butions for 8 GeV/c vr^ and p beams 



Fig. 2. Excitation energy (left frame) and residue mass probability (right frame) for 
several of the systems studied in this work, as indicated on figure. 



Fig. 3. Bottom: the probability for observing events with excitation energy greater 
than E* > 400 MeV relative to the probability for events with E* = 400 MeV. 
Systems are indicated on figure. Top: average mass loss AA in the fast cascade as 
a function of excitation energy. Systems are defined in bottom frame. 



Fig. 4. Probability distributions for data in Fig. 2, plotted as a function of E*/A of 
the residue. Systems are defined on figure. 



Fig. 5. Top: average ratio of observed and geometry corrected IMFs per residue 
nucleon as a function of E*/A; symbols are defined in bottom frame. Bottom: 
power-law parameters r from fits to the charge distributions as a function of E*/A 
of the residue. Systems are defined on figure. 

Table 1 

Ratio of the integrated events beyond multifragmentation threshold to total events 

with E* > 400 MeV {E*/A > 2 MeV). 
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